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Abstract 

A model for the symmetric coupling of two self-oscillators is presented. The non- 
linearities cause the system to vibrate in two modes of different symmetries. The 
transition between these two regimes of oscillation can occur by two different scenar- 
ios. This might model the release of vortices behind circular cylinders with a possible 
transition from a symmetric to an antisymmetric Benard-von Karman vortex street. 
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Nowadays the understanding of self-oscillators is fairly complete thanks to the 
theory of bifurcation and of normal forms. A familiar model for this is the van 
der Pol system Q that displays a wide range of behavior, from weakly nonlinear to 
strongly nonlinear relaxation oscillations, making it a good model for many practical 
situations. However there are physical situations characterized by spontaneous self- 
oscillations with certain basic features that are absent from the "generic" van der 
Pol system. Take for instance the Benard-von Karman vortex street in the wake of 
a cylinder. Its phenomenology is approximately as follows || |3|, [| : the velocity field 
remains two-dimensional for Reynolds number (Re) less than 160 (creeping flow for 
Re < 4 ; recirculation zone with two steady symmetric eddies attached behind the 
cylinder for 4 < Re < 45 ; instability at Re ~ 45 at wich these eddies are released 
alternatively to form a double row of opposite sign vortices, the Benard-von Kar- 
man vortex street) and for Re > 160 three-dimensional and irregular fluctuations 
are superimposed on the dominant periodic vortex shedding. 

It is tempting to say that the periodic vortex shedding provides a classical example 
of the Poincare-Andronov bifurcation to a limit cycle. However, one fundamental 
ingredient would be missing if one insisted in describing these oscillations by the van 
der Pol equation: no equivalent of the symmetry of the system would be present 
in this mathematical description. That is, this mathematical picture would make 
no difference between a symmetric and an antisymmetric release of vortices, both 
would be fairly described by the same van der Pol equation, although they are clearly 
physically different. 

We propose here to implement the major symmetries of the Benard-von Karman 
oscillations by assuming that they result from the symmetric coupling of two iden- 
tical oscillators, each one responsible of the periodic release of vortices on one side 
of the cylinder. The interest of this approach is that it shows two possible stable os- 
cillating states: one symmetric, one antisymmetric, depending on the value of some 
coupling parameter. By varying continously the coupling, it is possible to monitor 
the transition between these two regimes, something that is beyond an approach 
using a single van der Pol equation. 
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A dynamical model for representing these properties is presented and some in- 
formations about the transitions between the two oscillation modes are obtained. 
We make no attempt to relate our model to the fluid mechanical equations. The 
symmetry properties of the system are used as a basic ingredient, as well as the fact 
that it operates in a stable way in an oscillating mode. 

Let us assume that there is an oscillator (the 'vortex' emitter) on each side 
(sidei,side 2 ) of the cylinder H and let (xi,x 2 ) be their amplitudes of oscillation. If 
X2(t)) is a possible dynamics, (x2(t),xi(t)) is also realizable by symmetry. If 
a vortex is emitted from sidei whenever Xi (t) reaches a maximum then the symmet- 
rical and antisymmetrical vortex streets (Fig. 1) appear as two oscillation modes, 
one with the two oscillators in phase and the other one with the two oscillators out 
of phase. The symmetric mode Gi verifies xi(t) = X2{t) and the antisymmetric one 
2 verifies Xi(t) = x 2 (t + T/2) with T the period. 

The simplest model representing these properties is a system of two coupled har- 
monic oscillators with a small coupling f3: 

X\ + X\ + [3x2 = 

x 2 + x 2 + (3x 1 = (1) 

The normal modes verify: 0j + (1 — ( — l)*/3)Oj = (i = 1, 2), with 0i = X\ + X2 
and 9 2 = Xi — x 2 , each one with frequencies lo\ — 1 + (3 and u) 2 = 1 — (3. But this 
model is Hamiltonian and it is not useful to describe self-oscillations (i.e. oscillations 
resulting from balance between energy input and dissipation). It does not present 
any Poincare-Andronov bifurcation, although it was shown experimentally that the 
vortex shedding behind a cylinder results from that type of bifurcation To 
remedy this, we can introduce, as in the van der Pol system, a nontrivial damping 
term xf(x) in equations (1). The van der Pol oscillator of equation x— e(l— x 2 )x+x = 
0, is probably the simplest example of a system with one stable limit cycle: a fixed 
point at the origin and an unstable closed orbit for e < and an attractive cycle for 
e > 0. A natural extension of the van der Pol system to the representation of two 



3 



symmetric coupled oscillators is: 



Xi — e(l — (x\ + x\))±i + x 1 + f3x 2 = 
x 2 - e(l ~ (x\ + x\))x 2 + x 2 + (3xi = 



(2) 



The small displacements near x\ = x 2 = x\ = x 2 = are damped to zero when 
e < and give sustained oscillations when e > 0. The birth of a stable limit cycle 
is then governed by the parameter e and our interest rests in this regime (e > 0). 
Other properties of (2): 

(i) If /3 — 0, equations (2) present an 0(2) symmetry. If the complex variable 
z = Xi + ix 2 is defined, system (2) becomes: z — e(l — \z\ 2 )z + z = 0, which has the 
symmetries: z — > e^z and z — > z. A stable solution is z = e^r(t) with constant 
and r (i) solution of: f — e(l — r 2 )r + r = 0. Its representation point is a straight line 
through the origin in the (xi,x 2 ) plane at constant angle 0, r(t) oscillating along 
this line. The periodic solution z = e ±lt is unstable. 

(ii) If (5 ^ the phase symmetry is destroyed, although a Z{2) symmetry z — > ±iz 
remains. The equations (2) become: 



Now there are two oscillating solutions: the symmetric mode ©i = (x\ = x 2 ) 
given by Z\ = e* 7r / 4 r 1 (t) and the antisymmetric mode 2 = (xi = —x 2 ) given by 



z 2 = e m / 4 r 2 (t), where ri(t) and r 2 (t) verify: fj — e(l — rf)r,i + (1 — (— l) l /3)rj = 



This system presents stable oscillations when < 1 and diverges to infinity 
when |/3 1 > 1 (except for initial conditions r — f — 0). One finds a parameter 
value f3 e > such that: if < (3 < (3 e then Oi is stable and 2 unstable and, if 
— f3 € < j3 < 0, Oi is unstable and 2 stable. Also, if f3 € < \/3\ < 1 the two modes 
0i and 2 are linearly stable. Let us remark that this simple model brings all the 
information we are looking for. The range of parameters e > would modelize the 
situations of stable limit cycle oscillations observed experimentally for for Re > Re c , 
where Re c is the Reynolds number at the onset of vortex shedding. The parameter 



z-e(l 



z\ 2 )z + z + i(3z = 



(3) 



(i=l,2). 
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(3 would represent for instance the aspect ratio in the experiments of Le Gal and 
collaborators ||. 

In model equation (2) the transition from 0! to 02 stable oscillation occurs at 
(3 = 0. In this case the coupling is lost and the system becomes degenerate at 
transition (the phase difference between the two oscillators is arbritary) presenting 
an infinity of stable oscillating states. In order to remove this degeneracy, we need 
to have more than one coupling parameter. This means that the dimension of 
parameter space for a transition between two modes of oscillation should be greater 
than one: the unfolding of this transition should be controlled by two parameters at 
least. Thus, a more general and 'robust' scheme of transition from mode 0i to mode 
02, and viceversa, is achieved by introducing another phase symmetry breaking term 
(proportional to 7) in the dissipative force of equations (2): 

xi — e(l — x\ — (1 + i)x^))x\ + X\ + (3x 2 = 

x 2 - e(l - x\ - (1 + j)x{))x 2 + x 2 + fix! = (4) 

where 7 and (3 are the coupling constants. Symmetries (x\,x 2 ) <-> (x 2 ,xx) and 
(xi,x 2 ,/3) <-> (x 2 , —Xi, —(3) are preserved. 

We have numerically studied the solutions of this system in two different regimes 
and found the following results for 7 and (3 near zero : 

(a) When \(3\ <C I7I there are four oscillatory states: the pure symmetric mode 
X = (x\ = x 2 ), the pure antisymmetric mode ©2 = {xi = —x 2 ) and two new mixed 
modes ©12 = (^1,^2) and ©21 = (x 2 ,xx) intermediate between ©i and ©2- If 7 > 
the mixed modes are stable and the pure modes unstable. If 7 < the mixed modes 
are unstable and the pure modes stable. 

(b) When \(3\ 3> I7I equations (4) tend to equations (2) (the perturbation introduced 
by 7 can be neglected in front of (3), the mixed modes ©12 and ©21 collide and 
disappear, and the pure mode ©i and @ 2 remain. 

Let us explain in more detail the two different scenarios (Fig. 2) that can be 
found for the transition between the pure modes ©i and ©2 when 7 is fixed and (3 
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is varied (e is kept constant and of order 1, but the results are not sensitive to its 

specific value): 

Scenario I ; 7 < (Fig. 2a): 

(11) (3 < — c(e)|7| (c(e) positive constant, depending on e and of order 1 for e of order 
1): 0i is unstable and B2 stable. No mixed modes. 

(1 2 ) — c(e)|7| < P < c(e)|7|: the two unstable mixed modes 612 and 6 2 i grow from 
0i for (3 = —c(e)\j\. In this regime the two pure modes are stable. Depending 
on initial conditions the system oscillates in the symmetric or in the antisymmetric 
mode. When (3 — > c(e)|7| the two mixed modes approach 2 , collide with it for 
(3 = c(e) I'yl making 2 linearly unstable. It transfers the stability from 2 to 0i. 

(13) (3 > c(e)|7|: ©i is stable and 2 unstable. No mixed modes. 
Summarizing: there is a range of parameters J 2 where ©i and 2 are both stable, 
and each mode ( ©i or 2 ) loses its stability by a supercritical bifurcation on the 
edges of J 2 . 

Scenario II ; 7 > (Fig. 2b): 

(Hi) (3 < — c ( e )|7l : ©i is unstable and 2 stable. No mixed modes. 

(11 2 ) — c(e)|7| < (3 < c(e)|7|: the two mixed modes ©i 2 and 2 i bifurcate from 2 
for (3 = — c(e)|7|. These are stable (which makes the difference with scenario I). 
In this II2 regime the two pure modes are unstable and the system will decay in 
one of the two mixed modes according to the initial conditions. When (3 —>■ c(e) J'y] 
the mixed modes approach ©i and collide with it for (3 = c(e)|7|. It transfers the 
stability from the mixed modes to ©i. 

(113) (3 > c(e)|7|: ©i is stable and @ 2 unstable. No mixed modes. 

Sumarizing, there is a range of parameters 7/ 2 where the two mixed modes are 
stable, and collide with ©1 or 2 on the edge of II2 to exchange stability. 

A derivation of the dynamics of equations (4) can be obtained in the formalism of 
a slow phase dynamics 0. (A different calculation can be found in reference || where 
two nonlinear oscillators with diffusive coupling, not the one we consider, are studied 
in the vicinity of a Hopf (Poincare-Andronov) bifurcation). When (3 = 7 = 0, the 
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set of equations (4) present phase symmetry (0) and temporal translation symmetry 
(ijj), and when (3 or 7 are different from zero the phase symmetry is broken. For 
small P or 7 (of the same order of magnitude) the general solution of equations (4) 
can be written: 

x l = r o (t + ip) cos + 5x\ 
X2 = r (t + -0) sin + 8x 2 

where r (t) is the periodic non zero solution of the van der Pol equation: f — e(l — 
rl)f + r = 0, ip and follow a slow dynamics (^, ^ <C and 5x2, "0, are 
small and of order (7, (5). Linearising equations (4) to order (7, /3) we obtain a set 
of coupled equations to be solved for 8x\ and 5x 2 . These equations are written in 
matrix notation to make their structure more transparent: 

C ( f 1 ) = ( f t?\ Sin t 9( f°\™t ) ( t V^+jMro,^)] f ""I 1 (5) 
\Sx 2 J \ -/(r o )cos0 5f(r o )sm0 ) \ip ) L 2 /J \ cos0 y 

where 

£- ' ) + W)("f' ^) 

^ = 9„-e(l-r 2 )9 t + l 

/(r ) = 2r -e(l-r 2 )r 

fl^o) = -[2r - e(l - rjj)r ] 

h(r , 0) = er r sin(20) 

The relevant solution of equation (5) is made of periodic functions of time, with the 
same period T as r (t). This excludes functions with a secular growth and leads 
to a solvability condition that will ultimately become an equation of evolution for 
0(t). To write this solvability condition, one needs to define first an inner product 
of functions of time with period T as: < 6\a >= J T (^iO"i + 02(T2)dt, (0 = (^1,^2) 
is written as a two component vector). One notices now that the linear operator C 
has a nonempty kernel: 

^ / sin0 \ _ . / cos0 \ 

£cj = =>- w a = r , w b = r . , 

y — cos J y sin y 
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Because of this non empty kernel, the equation (5) has no solution in general that 
is periodic with period T. To have such a solution, the right hand side (p of this 
equation must be orthogonal to the kernel of the adjoint operator £ + , made of two 
functions, Xi (i — a,b), oi t that are solution of the formal equation C + x = 0. The 
solvability condition is then that the two inner products < Xi\<P > (i — a,b) are 
zero. The operator C + can be written explicitely as: 

^(f;)*«(r t ;j 

where JF+ = d tt + e<9 t (l — r%) + 1. Since the two left vectors Xa,b, once multiplied 
with the inner product, < | >, with the left side of equation (5) give zero, the same 
product with the right side of (5) should give zero as well. This gives two coupled 
equations for and tp: 

(t t)(lH(::M:) <«» 

where ha, h i2 , rrii, rii (i = a,b) are functions of e and (j) after the time integration 
coming from the scalar product: 

[T 

fki = / /(r o )[sin0 xn ~ cos0 Xi2\dt 
Jo 

h i2 = / ^(r o )[cos0 xn + sin0 xa]dt 
Jo 

f T 

rrii = r {sm(j) xn + coscj) Xi2]dt 
Jo 

1 f T 

rii = - / r h(r , 0) [sin xn + cos xa]dt 

2 Jo 

The 0-dependance of the vectors in the kernel of C + can be factored out by noticing 
that these vectors can have the following 0-dependance: 

From this the two functions h a ^(t) are the non trivial (=non zero) solutions of period 
T of the two linear homegeneous equations: 

F+[h a (t)] = [d tt + ed t (l-r 2 ) + l]h a (t) = 
[^ t + 2e(l-r 2 )d t ][h b (t)} = [d tt + e(l-r 2 )d t + l}h h (t) = 

8 



Then the equation for <fi is simplified to: 

&a(e)0 = / a (e)cos(20)/? + s a (e)sin(40) 7 



(7) 



where k a , l a and s a are functions of e only that are proportional to various scalar 
product of functions on the equation (6) with h a (t). Thus, k a (e) — — Jo f{To)h a {t)dt 
! la(e) = J T r h atb (t)dt ; and s a (e) = f J T r%r h aib (t)dt. 

The equation (7) presents, as j3 and 7 vary, the same bifurcations as the one found 
numerically for the original set of equations (4). To show this, let us define (3' = 
fij^fe) anc ^ 7* = ^f(e)' ^ na ^ wm ^ e cons idered now as the bifurcation parameters 
(the quantities and |44 are constants of order 1 at a fixed finite value of e, and 
so can be get rid off by scaling). The fixed points of (7) are roots (in 0) of: 

/3'cos(20) + ^-sin(40) = (8) 

or of cos(20) = or f3' + 7'sin(20) = 0. 

If \P'\ > |y|, this corresponds to scenario 7i j3 and 1/1,3. The only steady states are 
at the zeroes of cos(20), that are at = | and = — with one stable and the 
other unstable, depending of the sign of (3' (and consequently of (3) in agreement 
with what was found numerically. If \(3'\ < |7'|, they are two more fixed points, that 
are |sin~ 1 ( — ^) and | — | sin _1 (— They correspond to the mixed modes and, 
as j3' goes for instance from —7' to 7' (if 7' > 0), one finds the same bifurcation 
structure as found for the original equations (4), as explained formerly under the 
heading "Scenario II" (Fig. 2). 

In this rapid communication, we have presented a simple model for systems made 
of two symmetric coupled self-oscillators || . This might be a theory for one of the 
most studied instabilities in real fluid mechanics, the periodic release of vortices in 
the wake of cylinders, a phenomenon studied experimental and theoretically long 



ago by H. Benard and T. von Karman |L0], |TTJ and his collaborators. The connection 



of the present work with the Benard-von Karman phenomenon could be as follows. 
Our idea is that the wake is created by two symmetrically coupled self-oscillators, 
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one on each side of the cylinder. We have shown that, depending on the coupling, 
these two systems may either oscillate in phase or out of phase (as in the Benard- 
von Karman wake in a normal viscous fluid). Moreover, the transition from one of 
these two states to the other is realized by two different scenarios depending of the 
parameters. This might describe recent experiments by Le Gal and collaborators 
M, who observe this transition when the flow around the cylinder is more and more 
constrained by plates perpendicular to the axis of this cylinder. 
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Figure Captions 



Fig 1. A graphical representation of the two oscillation modes: (a) the symmet- 
ric vortex street Gi and (b) the antisymmetric one G2. 

Fig 2. Nonlinear transition between the two oscillation modes (61,62): (a) in 
scenario I (7 < 0) the two intermediary mixed modes, 0i2 and G21, are unstable and 
(b) in scenario II (7 > 0) these mixed modes are stable, (c) Another representation 
of scenarios I and II (inspired from figure 3 in reference ||). 
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